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Abstract 

Herein  the  electromagnetic  scattering  is  determined  using  the  Finite  Element  Method.  In  particular 
the  radar  cross  section  of  the  scatterer  is  estimated.  Comparison  between  the  Finite  Element  Method 
(FEM),  Method  of  Moments  (MoM)  and  the  method  of  Physical  Optics  (PO)  is  made  in  one  numerical 
example.  In  another  numerical  example,  the  convergence  rate  is  compared  using  pure  /i-refinement, 
7>enrichment  and  with  a  priori  /ip-adaptive  refinements.  The  gain  using  the  /ip-adaptive  approach  is 
significant  since  an  exponential  convergence  may  be  obtained. 

1  Introduction 

Accurate  calculations  of  the  electromagnetic  scattering  from  objects  is  of  great  importance  in  the  design 
of  low  observable,  so-called,  stealth  vehicles.  External  scattering  from  electrically  large  convex  objects 
is  usually  estimated  by  means  of  high  frequency  methods,  such  as  the  method  of  physical  optics.  So 
called  numerically  exact  methods,  such  as  method  of  moments,  finite  difference  and/or  the  finite  element 
method  are  preferable  whenever  it  is  possible. 

Here  the  electromagnetic  scattering,  expressed  in  terms  of  radar  cross  section  (RCS),  is  calculated 
using  an  /ip-adaptive  finite  element  method.  The  /ip-adaptive  method,  where  the  elements  have  locally 
variable  order  of  approximation  p  and  element  size  h,  is  very  suitable  for  modeling  singularities  such 
as  those  resulting  from  irregular  geometries.  Also,  the  numerical  dispersion  error,  e.g.  in  resonating 
cavities,  may  be  reduced  using  the  /ip-method.  The  ability  to  model  singularities,  is  one  of  the  largest 
advantages  for  using  the  /ip-method,  since  the  maximal  convergence  may  be  retained.  The  /ip-adaptive 
FE  code  used  here,  called  EM3D,  is  a  derivative  of  a  general  package  for  /ip-adaptive  FE  discretizations 
called  3Dhp90,  see  [1].  EM3D  has  extensively  been  applied  in  cavity  problems,  see  e.g.  [2];  herein  focus 
will  be  made  solely  on  the  scattering  from  a  compact  body  in  the  free  space. 

2  Problem  statement 

Consider  the  following  formulation  of  the  time-harmonic  Maxwell’s  equations  in  terms  of  the  scattered 
electric  field  E  defined  as  E  :=  Etot  —  Emc,  where  Etot  is  the  total  electric  field  and  Einc  is  the  incident 
electric  field.  Let  9  be  a  open,  bounded  domain  UcK3,  with  boundary  T  surrounded  by  the  free  space 
:=  M3\P  .  Let  also  the  exterior  domain  S7+  be  truncated  by  introducing  an  artificial,  closed  boundary 
To,  and  denote  the  part  of  the  reduced  exterior  by  hip"  and  the  enclosed  remaining  part  of  by  Ho,  he. 
fl+  :=  Ho  U  Pq".  Then  the  problem  is  to  find  the  scattered  electric  field  E(x),  *60u  Po>  such  that 

V  x  xeJ-  k20erE  =  -jk0Z0Jirnp  (1) 

is  satisfied.  In  the  above  er  :=  er(l  —  j—)  and  /ir ,  er  and  a  are  relative  permittivity,  relative  permeability 
and  conductivity  respectively;  and  co  is  the  circular  frequency.  Also,  in  the  above  equation  ko  :=  u)y/e opo 
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is  the  wavenumber  in  free  space,  and  Zq  =  yVoAo  is  the  intrinsic  impedance  of  free  space.  In  Eq.  (1) 
jimp  -g  impressed  current  defined  as 

Jimp  :=  ju>e0Einc,  Einc(x)  :=  E0  p  e~j  k° d  x .  (2) 

In  the  above  Eq  is  the  amplitude,  p  is  a  unit  vector  identifying  the  polarization  of  the  incident  field,  and 
d  is  the  direction  of  propagation. 

At  a  perfect  electrical  conductor  (PEC)  surface  T  the  non-homogeneous  Dirichlet  boundary  condition 

nx  E  =  -nx  Einc ,  (3) 


is  imposed;  n  =  n(x)  is  the  unit  vector  normal  to  T  at  a  point  x.  At  material  interfaces  the  following 
continuity  conditions  are  applied 


nx[Ej  =  -nx  [. Eincj  , 


n  x 

—V  x  E 

=  —n  x 

'  1  | 

Mr 

.Mr  J 

V  x  Einc, 


(4) 


where  [•]  denotes  the  jump  in  the  argument  across  the  interface. 

Energy  conservation  require  the  resultant  scattered  field  E  to  decline  according  to  Silver-Miiller 
radiation  condition 

lim  r  [V  x  E  +  jko  er  x  E]  =  0,  (5) 

r— >oo 

as  r  — >  oo,  where  r  is  the  distance  between  the  center  of  the  object  and  the  actual  observation  point. 
Infinite  Elements  (IE)  are  used  here  (in  the  domain  11))")  to  implement  Eq.  (5). 

The  variational  formulation  of  the  above  strong  formulation  is  obtained  by  multiplying  Eq.  (1)  by 
a  vector  test  function  F  £  W*,  integrating  over  the  domain  fi,  integrating  by  parts  and  using  Eq.  (4). 
The  corresponding  weak  formulation  can  be  stated  as  follows.  Find  E  £  W  such  that  B(E ,  F)  =  L(F), 
where  the  sesquilinear  and  antilinear  functionals  B  and  L  are  defined  as 


B(E,  F)  = 

y={n,n0,n+} 


[  —(V  xE)-(V  xF)dV  -  fcg  [  erE-  F  dV  \  , 

!V  l  l  'r  JV 


L(F )  =  -jk0Z0  /  jimP  '  F  dV~  /  n  x  (V  x  Einc) 

y={n,n0,n+}  v  r 


1 

Mr 


FdS, 


(6) 

(7) 


VF  £  W*  satisfying  the  homogeneous  equivalents  of  Eqs.  (3)  and  (4).  The  bar  above  the  test  function, 
F,  denotes  the  complex  conjugate.  Here  W  is  the  space  of  trial  functions  defined  as 

W  :=  {F  £  L :  V  x  F  £  L?-1}er  x  (V  x  F)  -  jk0E  £  L2},  (8) 

and  W*  is  the  space  of  test  functions  defined  as 

W*  :=  {F  £  L2r  :  V  x  F  £  L2,  er  x  (V  x  F)  -  jk0E  £  L2},  (9) 

where 

L2±1  :=  {F  :  Fr  £  L2(H  U  U  Q^),rJzlFt  £  T2(H  U  U  H+)},  (10) 

with  Fr  and  Ft  being  the  radial  and  transversal  (tangential  to  a  sphere  centered  at  the  origin)  components 
of  the  field  F.  Note  that  a  necessary  condition  for  a  FE-discretization  of  the  electric  field  E  in  order  to 


stay  in  the  space  W  is  that  only  the  tangential  component  of  E  must  be  continuous  across  inter-element 
boundaries,  i.e.  W  £  Hcur\(Q) .  The  radar  cross  section  a  is  defined  as 

\Eoa(x)  ■  pobs | 2 


<7  :=  lim  4irr 

r— >oo 


(if) 

\Einc{x)\\2  V  ' 

where  E^x)  is  the  scattered  electric  far-field  ,and  pobs  is  the  observer  polarization.  The  scattered 
electric  far-field  may  be  obtained  from  the  near-field  by  the  following  transformation 

k0 


Eoo(x)  — 


4vr||a;| 


[x  x  M  +  Z0x  x  {x  x  J)e-jkor'*}dS{r ), 


(12) 


where  Tr  is  a  surface  enclosing  the  scatter,  x  :=  £c/|a;|,  r  is  the  current  point  on  the  integration  surface 
rr.  The  vector  fields  M  and  J  denote  surface  magnetic  and  electric  currents  defined  as 


M  :=  E  x  n,  J  :=  n  x  H , 


where  H  :=  x  E  is  the  magnetic  field. 


(13) 


3  Numerical  examples 

In  this  section  two  numerical  examples  will  be  studied. 
The  first  consists  of  an  open  pipe  where  the  electromag¬ 
netic  back-scatter  is  computed  using  FEM  (EM3D),  MoM 
(M-RANDOLPH,  which  is  a  apart  of  the  suite  of  programs 
within  GEMS,  see  [3])  and  PO  (FOPOL).  The  second  ex¬ 
ample  is  a  convergence  study  of  the  computed  scatter  from 
an  infinitesimal  thin  PEC  disc. 

3.1  Open  pipe 


Here  the  electromagnetic  scattering  on  an  open  pipe  will 

be  calculated.  The  dimensions  of  the  pipe  is  36.0  x  4.0  x  Figure  1:  RCS  at  1  GHz  (HH-polarization) . 
0.065  inches  (Length  x  Diameter  x  Thickness).  The  pipe 

is  considered  as  a  perfect  electrical  conductor.  The  mono-static  RCS  will  be  estimated  in  an  azimuth 
sweep  from  0  to  180  degrees,  where  the  x-axis  is  parallel  with  the  axis  of  the  pipe.  The  results  can  be 
seen  in  Figs.  1  to  3  for  three  different  frequencies.  At  an  angle  perpendicular  to  the  manifold  of  the  pipe 
all  three  methods  agree;  however,  when  the  direction  of  the  impinging  wave  approaches  the  direction 
of  the  x-axis  the  error  using  PO  is  significant.  At  lower  frequencies  MoM  agree  well  with  FEM,  but  at 
higher  frequencies,  where  more  energy  is  able  to  enter  the  pipe,  there  are  some  discrepancies  between 
the  methods. 


3.2  PEC  disc 

In  this  section  the  convergence  properties  of  separate  and  combined  polynomial  enrichment  and  mesh 
refinement  methods  are  studied. 

Pure  /r-refinement  and  pure  p-enrichment  convergence  rates  are  compared  with  /ip-refinement.  The 
object  under  investigation  is  a  PEC  disc  with  infinitesimal  thickness  and  radius  a  =  1.  For  small 
wavenumbers  an  asymptotic  series  solution  is  available.  The  wave  number  considered  here  is  ka  =  0.5. 

In  the  case  with  pure  h  and  p  refinement /enrichment  algebraic  convergences  are  obtained.  As  can 
be  seen  in  Fig.  4,  the  increase  in  slope  of  the  convergence  rate  from  using  pure  p-enrichment  instead 
of  pure  /i-refinement  is  approximately  1.7.  In  the  case  of  /ip-refinement  we  obtain  almost  exponential 
convergence  and  far  better  accuracy  than  with  pure  ^.-refinement  and  p-enrichment. 


Note  that  the  property  of  interest  here  is  the  Zo-norm 
of  the  RCS.  The  reference  solution,  denoted  by  aex,  is 
obtained  using  a  low  frequency  analytic  series  expansion, 
see  [4],  and  aje  is  the  calculated  FEM  solution.  Also  note 
that  the  initial  mesh  has  to  be  designed  so  that  the  wave 
is  resolved,  i.e.  it  has  to  fulfill  the  sampling  theorem.  The 
a  priori  hp- mesh  was  constructed  using  three  h  and  one  p 
refinement  in  each  sequence,  i.e.  the  last  /ip-refinement, 
in  Fig.  4,  consists  of  three  sequences  with  three  h  and  one 
p  in  each  sequence. 


Figure  3:  RCS  at  4  GHz  (HH-polarization). 

4  Conclusion 


Figure  2:  RCS  at  3  GHz  (VV-polarization). 


In  this  paper  the  electromagnetic  scattering  was  mainly  Figure  4:  Scattering  from  a  PEC  disc,  ka  =  0.5. 
determined  using  an  /ip-adaptive  version  of  the  FE-  Far-held  convergence  of  the  ^-norrn  of  the  RCS. 
Method.  When  considering  a  problem  with  existing 

singularities  it  is  shown  herein,  using  an  infinitesimal  thick  disc,  that  the  /ip-approach  is  superior  to  tra¬ 
ditional  approaches  such  as  pure  /i- refinement  and  p-enrichment.  The  convergence  rate  of  the  calculated 
RCS  is  almost  exponential  for  the  /ip-approach. 
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